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Unidirectional transports have been observed in two-dimensional systems, however, so far they have not been 
experimentally observed in three-dimensional bulk materials. In this theoretical work we show that the recently 
discovered Weyl materials provide a platform for unidirectional transports inside bulk materials. With high 
experimental feasibility, a complex Dirac mass can be generated and manipulated in the photonic Weyl crystals, 
creating unidirectionally propagating modes observable in transmission experiments. Possible realization in 
(electronic) Weyl semimetals is also studied. We show in a lattice model that, with a short-range interaction, the 
desired form of the Dirac mass can be spontaneously generated in a first-order transition. 

PACS numbers: 73.43.-f,78.67.Pt 


Introduction. Recently, novel materials described by the 
Weyl equation are among the research focuses in physics^-. 
For instance, the phenomenon of chiral anomaly^ has been 
extensively investigated in Weyl semimetals^i^. The recent 
experimental discoveries of the photonic Weyl crystals^— 
and the TaAs class Weyl semimetals^T^? have generated in¬ 
tense interests in Weyl materials. 

A Dirac mass can couple two Weyl points together to form 
a massive Dirac fermion (FiglT^)21. The Dirac mass is gener¬ 
ally complex-valued, and it is known in continuous field the¬ 
ory that the vortex lines of mass carry unidirectionally prop¬ 
agating modes^^. Since the unidirectional modes are com¬ 
pletely immune to backscattering, they are capable of support¬ 
ing dissipationless transport. Such unidirectional modes were 
found at the edge of quantum Hall systems and their photonic 
analogues^^, however, so far they have not been experimen¬ 
tally observed inside three-dimensional (3D) bulk materials^^. 
Therefore, it is highly desirable to seek experimental realiza¬ 
tion of the vortex line of Dirac mass, which is nevertheless 
challenging. In particular, designing a material with the phase 
of mass continuously tunable from 0 to 2n is a difficult tasfei^. 

The results of this paper are threefold. First, we study the 
unidirectional modes in 3D lattice Weyl materials, and con¬ 
struct a discrete version of the vortex line of mass. This study 
is complementary to the continuous theories mentioned pre¬ 
viously. New features absent in the continuous theories are 
analyzed. 

Second, we propose to realize the vortex line of mass and 
the associated unidirectional modes in photonic Weyl crys¬ 
tals. This can be achieved by a simple modulation of the 
photonic Weyl crystals, which naturally generates a complex 
Dirac mass, whose phase can be manipulated by lattice trans¬ 
lations. This proposal is most feasible after the recent ex¬ 
perimental realizationSi of the gyroid photonic Weyl crystals. 
This is our central proposal. 

Third, we point out that the spontaneous mass generation 
in Weyl semimetals can be first-order due to the lattice effect. 
This first-order nature is absent in continuous approximations 
with a momentum cutoff (e.g. Ref. —). 

Lattice model and Dirac mass engineering. Let us be¬ 
gin with a simplest lattice model, which captures the key 
low-energy physics of Weyl semimetals and photonic Weyl 



FIG. 1: (a) Coupling two Weyl points by a Dirac mass term, (b) A 
“mass clock” showing the four values of the complex Dirac mass in 
our lattice model. 


materials with a single pair of Weyl points. The lattice 
Hamiltonian^ is given as 

Ho{k) = 2tx sin kxCTx + 2ty sin kyCTy + [2l(cos kz - cos ko) 
+B{2-coskx-cosky)]o-z. (1) 

For simplicity, the lattice constants a has been scaled to 1, 
therefore k, has periodicity 2nla - 2n. We take the inter¬ 
pretation that each unit cell contains a single site, and the 
Pauli matrix cr, refers to the spin (or pseudospin) degree of 
freedom. This lattice model possesses two Weyl points at 
K± = (0,0, +ko), with difference Q = K_ - K+ = (0,0, Q) = 
(0,0,2ko). Near the Weyl points, we have the low-energy 
Hamiltonian 

Hwiod = VxCTxqx + VyO-y-qy + VzCTzTzqz, (2) 

where Vv 0 >) = '^Ix{y), Vz = 2tzSinko, the Pauli matrix is the 
chirality operator, namely, the eigenvalue of is +1 for Bloch 
states near K±, and k - K± s {qx,qy,qz) is the momentum 
measured relative to the Weyl points. Hereafter we shall take 
ko - njA, tx - ty - 0.5, tz - 1, B - 2.5, thus Q = (0,0,7r/2). 
The Dirac mass terms can be formally included as 

Hm — mcTzTy + m*(TzT-, 


(3) 
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FIG. 2: Lattice system with dislocation lines. Each small cube rep¬ 
resents a unit cell. Different colors represent different values of m, 
(namely, see text), (a) A vortex line of the complex 

Dirac mass (axion string), (b) A system with four vortex lines. Peri¬ 
odic boundary condition is taken in the x and y direction, while k- is 
kept as a good quantum number. The arrows indicate the propagating 
directions of the unidirectional modes. 


where T± s {Tx±iTy)l2, and both crjT+ andcrj.T- anticommute 
with the three terms in Hw- At low energy, the eigenvalues of 

the full Hamiltonian Hw+Hm are + ■^Yii=x,y,z 
a gap 2\m\. 

Now we would like to introduce a Dirac mass in the lat¬ 
tice model in Eq.([T]l. Because the momentum difference be¬ 
tween the two Weyl points is Q = (0,0, 7 r/ 2 ), a Dirac mass 
must break the translational symmetry, enlarging the size of 
unit cell by a factor of four in the z direction. The four adja¬ 
cent sites in each enlarged unit cell are sequentially denoted 
as A, B, C, D, namely, “A” for z - 4n, “B” for z = 4n H- 1, “C” 
for z = 4« H- 2, and “D” for z = 4« -H 3 (n being integers), where 
z refers to the z coordinate. To generate a Dirac mass, we add 
a translational-symmetry-breaking term into the Hamiltonian; 

H'^ (4) 

j 

where j runs through all sites, and nij - mA,mB,mc,niD for 
A, B, C, D site, respectively. The inclusion of cr^ in H' is sug¬ 
gested by the form of Hm- 

After enlargement of the unit cell, the lattice Hamiltoni¬ 
ans Hq and H' are replaced by eight-band lattice Hamiltonians 
Hq and H' (Supplemental Material), meanwhile, the Brillouin 
zone of the enlarged unit cell is [0, 2 ; 7 r] x [0, 2 ; 7 r] x [ 0 , 7 r/ 2 ]. 
We find that the energy eigenvalues of Hq have fourfold de¬ 
generacies at zero energy at k = (0, 0 , 7 r/ 4 ). For these four 
low-energy bands we obtain (Supplemental Material) the fol¬ 
lowing effective Hamiltonian of Hq + H': 

Hefiiq) - VxCTxqx + VyO-y.qy + v^a-^T^q^ + mcr{r+ + ina^T- 

(5) 

where {q^, qy, q^) = {kx, ky, k^) - (0,0, 7 r/ 4 ), Vx(y) = 2tx(y), = 



FIG. 3: The spectrum of the system illustrated in Fig[2}j. The system 
size is 28x28 unit cell in the xy plane, and nto = 0.5. Four dislocation 
modes inside the bulk energy gap can be clearly seen. (The tiny 
energy gap at 4k- = jt is due to finite-size effect.) Energy bands far 
from zero-energy are not shown here. 


2t^ sin( 7 r/ 4 ), A = (ma -i- + nic + mo)/4, and 

m — (ma + + e^'^mc + e^‘®m£))/4. ( 6 ) 

We shall study the A term later; at this stage we just omit it 
and focus on the Dirac mass m. 

Now we can construct a vortex line of the Dirac mass. 
We observe that for the four cases, (I) (m^, mg, me, mo) = 
(mi,m 2 ,m 3 ,m 4 ); (II) (m/t,mg,me,mo) = (m 4 ,mi,m 2 ,m 3 ); 
(III) (ntA, mg, me, nio) - (m 3 , m 4 , mi, m 2 ); (IV) 

(mA,mB,mc,mo) = (m 2 ,m 3 ,m 4 ,mi), the Dirac mass is given 
by mi = (mi-l-!m 2 —m 3 —!m 4 )/ 4 , mu = (m4+imi—m2—im-i)/4 — 
imi, mm = (m 3 -I- im^ — m\ — im2)/4 — —mi, and 
miv = (m 2 -I- /m 3 — m 4 — tmi)/4 = —imi, respectively. 
We can see that m is multiplied by a factor i during each move 
in I^II^III^IV^I. It is worth noting that configurations 
I, II, III, and IV are related by translations along the z 
directional. There is no difference in physical properties 
to distinguish them, unless they are adjacent to each other. 
Nevertheless, we can design a vortex line (an “axion string”), 
such that along a circle enclosing this dislocation, the path 
I—>11—»III—»IV—>I is followed (see Figl2^ ). This is a discrete 
version of the vortex line of Dirac mass (see Fig[T] 5 ). 

Hereafter we shall take (mi, m 2 , m 3 , m 4 ) = 

(mo, mo, -mo, -mo), which leads to A = 0. The Dirac 
mass for the four configurations is mi = exp()7r/4), 

mil = ™i, mm = —m^ and miv = —imi, respectively. 

In the continuous field theory (for instance, see Refi^), the 
vortex line of m carries a unidirectional mode. It is conceiv¬ 
able that our discrete version of vortex line of m also carries 
unidirectional mode. To be conclusive, we do a numerical cal¬ 
culation on a system with periodic boundary condition in the 
X and y directions (see FiglUt). In the calculation we have 
Fourier transformed Hq in Fq.([T]) into real space in the x and 
y directions. The energy spectrum (of the full Hamiltonian 
H = Hq 4- H') as a function of k^ is plotted in Figl3] from 
which we can see four unidirectional modes, two -Hz propa¬ 
gating and two -z propagating. 

To show that the unidirectional modes are indeed localized 
around the dislocation lines, we plot their density distributions 
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FIG. 4: The density distributions of the unidirectional modes prop¬ 
agating in the +z (Upper) and -z (Lower) directions, at momentum 
k, = 0.97;r/4. The parameters are the same as specified in Fig[3] 

in FiglU We only plot one mode for each direction of prop¬ 
agation (-I-Z and -z). Each mode is evenly distributed around 
two dislocation lines because of the tunneling effect (like in 
the double-well quantum mechanical system). 

Previously we have omitted the A term in Eq.(|5]l. With this 
term included, the energy eigenvalues are 

£'i,2,3,4(q) = ± ^jvlql + v]qj -I- (A + + |otP)2, (7) 

which become gapless when |A| > \m\. As long as |A| < \m\, 
the bulk material is gapped, and the topologically robust uni¬ 
directional modes persist. 

What if 27r/|Q| is incommensurate with the lattice constant? 
Let us take ko - n/4 + 6, with 6 being small. This changes 
A in Eq.(|5]l to A - 2 Lcos( 7 r /4 H- d) H- 2fj. cos( 7 r/ 4 ) » A H- 
(Supplemental Material). The previous condition |A| < \m\ is 
then replaced by |A -H v^d\ < \m\. Therefore, a small 6 does not 
qualitatively change our result. 

The above study is relevant to Weyl semimetals with 
time-reversal symmetry breaking, such as magnetic Weyl 
semimetals^^i^. In the next section we will study electron- 
electron interaction that may spontaneously generate the de¬ 
sired form of mass. More feasibly, we propose to realize the 
unidirectional modes in photonic Weyl materials, for which 
the Dirac mass can be readily generated and controlled. In the 
gyroid photonic crystals studied in Ref.— (Supplemental Ma¬ 
terial), one can tune Q to be close to (0,0, njl), and produce 
a modulation with a period of four unit cells (by making the 
A, B unit cells slightly different from the C, D unit cells). Erom 


symmetry consideration, this will generate the H' term at low 
energy, which leads to a Dirac mass. Then one can fabricate 
vortex lines (see Eigl2]l and observe unidirectional propaga¬ 
tion of light (or microwave) in transmission experiments. The 
discrete version of vortex line studied here is convenient to 
fabricate in experiments. 

Spontaneously generated mass. In this section we study an 
interacting electron model, whose free Hamiltonian is Hq - 
c^//o(k)ck, where Ck is the electron operator, and 

//o(k) = 2 l sin k^crx + 2fj, sin kyCTy + [2l(cos — cos ko) 
+B{cos kx - cos kj,)^]crj. (8) 

We shall take t^ - ty - 0.2, L = = 1.48, and 

k() - 7 r/ 4 . The Brillouin zone contains four Weyl points, 
K± = (0,0,+ko) and = (tt, tt,+ ko). The low energy elec¬ 
trons near K± can be described by Eq.(|2]i, while those near K'^. 
can be described by 

^^(*1) = -VxO-xqx - VyO-yqy + v^o-^T^q^, (9) 

In the continuous theory, it is well known that sponta¬ 
neous mass generation can be induced by a four-fermion 
interaction^, which suggests us to consider a Hubbard-like in¬ 
teraction H[ - -U 2;(cjcrjC;)^ in our model. It can be decou¬ 
pled as Yj Mjc^jcr^cj + UM^) at the mean-field level—, 

where we have introduced Mj = {c^cr^cj}. The mean-held 
Hamiltonian is therefore 

Hmp ^Ho- 2uY^Mjc)(r,Cj, ( 10 ) 

j 

Our mean-field calculation shall be quite standard, except that 
we shall enlarge the unit cell to x Ly x L, to accommodate 
possible symmetry-breaking phases ( LxXLyXL, = 1x1x4, 
2x2x4, 1x1x8, and 2 x 2 x 8 are calculated). Since 
Hmf is quadratic, it can be diagonalized in the k space, with 
eigenvalues e„(k), and the mean-field energy per unit cell is 
Emf = f ■(03 Ya fa(k) + U Yi=i’ where the summation 
Ya is over all negative energies eQ.(k) < 0 , i.e. the energies 
of the occupied Bloch states. The mean-field ground-state can 
be obtained by minimizing Emf with respect to M,, leading 
to standard integral equations^S for Mi (i - 1 ,2, ■ ■ ■, L^LyL,), 
which can be solved iteratively in numerical calculations^. 

For the 1x1x4 unit cell, the obtained values of 
(Ml,M 2 ,M 3 ,M 4 ) are plotted in Fig|5l All M,’s are zero 
for small U. A.t Uc ~ 1.12, there is a first-order tran¬ 
sition to a symmetry-breaking phase with M, 4^ 0. The 
same phase is reproduced in calculations on larger unit 
cells. For instance, we find that the values of (M,) are 
(Ml, M 2 , M 3 , M 4 , Ml, M 2 , M 3 , M 4 ) for the 1x1x8 unit cell, 
with the same values of Mi ,2,3,4 as obtained from the 1x1x4 
unit cell. A first-order transition is understandable here be¬ 
cause £’mf(M/) 4 £’mf(-M,), so that Emf contains odd- 
order terms when expanded in powers of M, ’s (similar to the 
Ginzburg-Landau theory of first-order nematic-to-isotropic 
transition in liquid crystals). 

Comparing Eg.lfTOb with Eq.®, we can see that the physics 
of the previous section applies here if we identify mj - 
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FIG. 5: The ground-state M;’s as functions of the interaction param¬ 
eter U . The Dirac mass is nonzero when U > Uc (see text). 



FIG. 6: Upper: The spectrum of the system with dislocation lines 
illustrated in Fig[2j), for the mean-field Hamiltonian Hmf at (7 = 
1.125. The system size is 28x28 in the xy plane. Lower: The density 
distribution of one of the unidirectional modes (at k, = 0.97;r/4). 


—2UMj\ in particular, the Dirac mass ism - -^(Ma + iMs - 
Me - iMj)). Dislocation lines can be constructed in similar 
way to the previous section. The spectrum and the density 
distribution of one of the eight unidirectional modes is shown 
in Figl6](Due to the presence of two copies of Dirac fermions, 
the number of unidirectional modes is twice that of the previ¬ 
ous model). 


Finally, we remark that dynamical mass generation is not 
the only approach to create a Dirac mass. In the Weyl 
semimetals based on layered structures^, we suggest to in¬ 
troduce an additional superlattice structure with period near 
27r/|Q| (say a modulation of the magnetization with period 
near 27r/|Q|), which will generate a Dirac mass. Then vortex 
lines of the Dirac mass can be constructed (similar to FiglU. 


Final remarks. We have shown how to create unidirectional 
modes in the bulk materials of electronic and photonic Weyl 
materials. Compared to the helical dislocation modes in weak 
topological insulators^^, our unidirectional modes are com¬ 
pletely immune to any kind of backscattering. In this respect 
our proposal will be of significant practical use. 


We have proposed the photonic Weyl materials as the most 
experimentally feasible systems, because the Dirac mass can 
be readily designed and controlled therein. The recent exper¬ 
imental observation of photonic Weyl points^S is most propi¬ 
tious to this proposal. Finally, we remark that our proposal 
may be generalized to other topological materials (e.g. topo¬ 
logical phononic or acoustic materials^—). 
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This can be achieved by spontaneous mass generatio n"i‘’4‘|S8- ^^ 
’s For conceptual discussions, see ReffiliS 

’s Vortex line of Dirac mass in the Weyl semimetal has been briefly 
discussed in Reffli, though no concrete model was studied there. 
” In the Dirac semimetalsSii, there is no such a relation between 
the Dirac mass and lattice translations. 

Analogous to Eq.(3.45) of Ref4^. 

In the cases of several solutions, we should take the one with low¬ 
est mean-field energy as the ground-state. 


Supplemental Material 


I. DIRAC HAMILTONIAN AND DIRAC MASS IN THE 
ENLARGED UNIT CELL 

In the paper, we have seen that the translational-symmetry- 
breaking term H' = enlarges the size of the unit 

cell by a factor of four in the z direction [Note: cr^(j) is the 
Pauli matrix <t^ defined at the j site]. In the enlarged unit cell, 
we have 2x4 = 8 bands. Now we derive the low-energy Dirac 
Hamiltonian and the Dirac mass. 

For the enlarged unit cell, the first Brillouin zone is [0,2;7r] x 
[0,2;7r] X [0,7r/2] (the lattice constant a of the original unit 
cell is taken to be 1). According to the standard band the¬ 
ory, within the description of the enlarged unit cell, the lattice 
Hamiltonian Hq becomes 


Hoik) = 


Hn 

Hn 

//3i 

//41 


Hn 

Hn 

Hn 

Hn 


Hn 

H 23 

Hn 

Hn 


Hl4^ 

H 24 

H 34 

H 44 


( 11 ) 


where //, ,’s are 2x2 matrices: 


Hi 1 = H 22 - Hn - H 44 - 2U sin + 2ty sin kyCTy 
-i-[-2fjCosko -H 8(2 - coskx - coskj,)]^^, 

Hn - H 21 = H 23 - Hn - Hn = H 42 - t^cr^, 

Hi4 = hI^ = exp(-4ik^)a-^, 

Hn = Hn = H24 = H42 = 0 . ( 12 ) 

The eight components of the state vectors have been ordered 
as (A t, A i,B '\,B i,C t, C J,), t and i referring to 


ko - 7r/4, U = ty = 0.5, B - 2.5. The only gapless 

point in the Brillouin zone is (0,0,7r/4) [Note that (0,0, -njA) 
is the same point because has periodicity 7r/2]. We can find 
that Hq have four zero eigenvalues at (0,0,7r/4), the corre¬ 
sponding eigenvectors are 


I'Ai) = 7 : 



f ^ 1 


f ® 1 


0 


1 


e-ml4 


0 

1 

0 

1 

e-ml4 

2 


103 

II 

0 


0 


e-ml2 


g-3iJr/4 


0 


. 0 , 


1 e-3,V4 j 


I<A3) = 7 : 



f ^ 1 


f ® 1 


0 


1 


eml4 


0 

1 

0 

1 

em/4 

2 

em/2 

103 

II 

0 


0 


em/2 


g^m/4 


0 


. 0 , 




(13) 


The existence of four-fold degeneracy at zero energy is under¬ 
standable because the two Weyl points in the original small 
unit cell are folded to the same point in the Brillouin zone 
of the enlarged unit cell. Here, \ifi) and liki) have cr^ = +l 
and -1, respectively, and similarly for 11 ^ 3 ) and 11 /^ 4 ). For later 
convenience, we also introduce a quantum number t^, which 
is 3-1 for and \ 4 i 2 ), -1 for and \fi 4 ). We shall not be 
concerned with the other four eigenvectors with energies of 
order of t^. 

Near the gapless point (0,0,7r/4), for the four low-energy 
bands, we can use the standard k ■ p perturbation theory^ to 
obtain a low-energy effective Hamiltonian i?o,eft, whose en¬ 
tries are 


(Ho,eff(tl))(40 ~ ^ VPp) 3" 

i=x,y,z ' 


(14) 


where iqx,qy,qz) = (kx,ky,kz) - (0,0,7r/4), a,(5 - 1,2,3,4, 
the partial derivatives are taken at (0,0,7r/4), and “■ • •” refers 
to higher order terms, which are omitted in the low-energy lin¬ 
ear approximation. These matrix elements can be calculated 
explicitly. For instance. 


dHo 

dkz 


0 

0 

0 

-AitzO-z 


AitzO-z 3 
0 
0 
0 


(15) 


therefore, we have 

,dHo 


(lAil 


dkz 


Itkl) = -<lA2l^^l<A2) = 

okz okz 




fiJ-f 

-<iA4l^|iA4) = 24sin(;r/4), 
ok. 


( 16 ) 
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while = 0 for other a,/l pairs. All other matrix 

elements can be calculated similarly, and the hnal result can 
be written concisely as 

Ho,es = + VyqyO-y + v^q^cr^T^, (17) 

where Vx(y) - ^txfy), - 2t^ sin(7r/4). This is a massless Dirac 
Hamiltonian. 

Next we add the H' - Yjj >njcr^{j) term. In the eight-band 
description of the enlarged unit cell, this mass term reads 


mAO-y, 

0 

0 

0 ) 

0 

'WfiCTj 

0 

0 

0 

0 

mcCTj 

0 

0 

0 

0 

ITIdCT^ ^ 


(18) 


Near the gapless point of Ho.eft, namely, (0,0,7r/4), we would 
like to obtain its effective form within the four low-energy 
bands, denoted by H'^g, whose matrix entries are 


(KsU = (19) 


By explicit calculations, it is not difficult to hnd that 


= -<t^2|/7V2) = <iA3|/7V3> 

= -<^/r4|HV4> = -^(mA +mB + mc + nio), (20) 

and 




- —(inA + + e'^nic + 

- ^{mA + irriB - me - imo), 

- -—{ma -h -I- e‘”mc + 

- ——(mA + ims — me — imo), 

= ««A2|HV4»*. ( 21 ) 


All other matrix elements of {il/a\H'\tljp) vanish. Therefore, at 
low-energy the Dirac mass term H' takes the form 




! A 0 
0 -A 
m* 0 
0 —m* 


m 0 ' 
0 —m 
A 0 
0 -A, 


( 22 ) 


where 


—(m^t -H iniB — me — imo) 

"^(mA + e'^ms + e^'^me + e^'^mo). 


(23) 


and 


A = 


-(mA +mB+me + mo). 


(24) 



9x/4- .I 

77r/4.I-! 

57r/4-. I- ' : 

3'7r/4 '.I-1 ; 

7r/4‘-1 j 

0 27r ^ 

ib) 


FIG. 7: The evolution of the phase 6 of the Dirac mass around 
the vortex line, for (a) an ideal vortex line of mass, and for (b) the 
discrete version of vortex line of mass in our paper. Here 0 is the 
angular variable in the cylindrical coordinates (z, r,(f)) ( r = 0 is the 
location of vortex line ). The rough profiles of the evolution of 9 are 
similar in (a) and (b), which justifies regarding the vortex line in our 
paper as a discrete analogue of the vortex string in continuous field 
theory. 


More concisely, we have 

H'jf = mo-{T+ -I- m*cr^T- + Actj. ( 25 ) 

The full low-energy Hamiltonian is 

^eft = ^O.eff + ^eff 

= (VxqxO-x + VyqyO-y + v,q^a-{T^) 

+{mcr^T+ + m*cr^T-) + Acr^. ( 26 ) 

All terms anticommute with each other except the last term 
AcTj.. The physical effect of a nonzero A is to shift the location 
of the Weyl points, which does not significantly modify our 
results as long as |A| < \m\ ( see the paper for discussions on 
the effects of Actj., and a derivation of the condition | A| < |m|). 

As discussed in the paper, the Dirac mass for the 
four cases (I) (mA, ms, me, mo) — (mj, m2, m3, m4); 
(II) (mA, ffiB, me, mo) — {1114, mi, m2, m3); (III) 

(mA,mB,me,mo) — (m3, m4, mi, m2); (IV) 

{mA,mB,me,mo) — (m2, m3, m4, mi), is given by 

mi = (mi H- im2 — m3 — /m4)/4, mn = (m4 -H imi — m2 — /m3)/4 = 
imi, mni = (m3 -I- im4 — mi — /m2)/4 = —mi, and 

miv = (m2 + im^ — m4 — imi )/4 = —/mi, respectively. 
Translating the system by one site along the z direction 
generates a phase factor = i for the Dirac mass, namely, 
m im. 

In the paper, we have described the construction of a 
vortex line of mass. For the case (mi,m2,m3,m4) = 

(mo, mo, -mo, -mo), the evolution of the phase of the Dirac 
mass around the vortex line is shown in Figl^t. 

The phenomenon that a translation generates a phase for the 
Dirac mass has an analog in the Su-Schrieffer-Heeger(SSH) 
model (see Ref 4 ^ for an introduction). The SSH model de¬ 
scribes a one-dimensional chain with hopping f,4+i = f 4- 
(-1 )'m/2 between the adjacent sites / and / 4-1, in other words, 
for two neighboring links, the hopping amplitude is f 4 - m /2 
and t - mjl, respectively. The low energy theory of the SSH 
model is an ( 14-1 )D Dirac fermion, the Dirac mass being equal 
to m. Translating the system by one lattice site generates a 
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phase factor n for the Dirac mass, namely, m —> -m, there¬ 
fore, there is no difference in the bulk properties to distinguish 
the m > 0 and the m <0 configurations. Despite of this, if we 
create a domain wall between these two configurations, there 
will be zero modes living on it. In two-spatial dimensions, a 
related model on the honeycomb lattice has been studied^. 


II. EFFECT OF INCOMMENSURATE ko AND Q 

Now we discuss the effect of shifting ko away from njA, say 
ko - njA + 6, 6 being small. The basic idea is not difficult to 
explain. Suppose that ko — njA + 6,6 being small, and the full 
Hamiltonian of this system is denoted as Hs + H' (namely. Ho 
is replaced by Hs). We split the Hamiltonian as 


Hs + H' ^Ho + H' + (Hs - Ho), (27) 

and treat the Hs - Ho term as perturbation. Let us first ignore 
the Hs - Ho term, then the problem is just the case we have 
studied, and the explanation of the dislocation line as discrete 
vortex line is valid. As we have shown, there are unidirec¬ 
tional modes in this case. Then we take the Hs - Hq term 
as a perturbation to this Hamiltonian. Since the Hamiltonian 
Hq + H' has already an energy gap ~ 2\m\ in the bulk sam¬ 
ples, the effects of perturbation is insignificant as long as it is 
not too big. In accordance with this insignificance of pertur¬ 
bation in the bulk sample, we do not expect the topologically 
robust unidirectional modes to suddenly disappear when the 
perturbation is small. Moreover, according to Callan and Har¬ 
vey, the existence of unidirectional modes is consistent with 
the anomaly inflow in the bulk material; since an insignif¬ 
icant perturbation in the bulk material should not eliminate 
the anomaly inflow, it should not eliminate the unidirectional 
modes either. 

Now we give a more quantitative criterion. In the enlarged 
unit cell, nonzero 6 amounts to replacing //,, in Ho by 


which become gapless when |A + v^6\ > \m\. As long as 
|A -H v^6\ < \m\, the bulk sample remains gapped, and the topo¬ 
logically robust unidirectional dislocation modes are expected 
to persist. 


III. TWO-BAND LOW-ENERGY MODEL FOR THE 
GYROID PHOTONIC WEYL CRYSTALS 


In the paper we have studied a simplified two-band model 
and show that vortex lines of Dirac mass carry unidirectional 
modes. Since our central proposal is to realize unidirectional 
modes in photonic Weyl materials, we would like to show that 
this model is indeed relevant to the experimentally feasible 
materials, the gyroid photonic crystals. 

Gyroid photonic Weyl materials have been proposed in 
Ref.^’ and observed in Refi22. The low energy physics of 
these photonic Weyl materials can be described by a three- 
band model (see the supplemental material of Ref.^*). Here, 
we would like to reduce the three-band model to a two-band 
model resembling the one studied in our paper. 

In the presence of a magnetic field B, the low-energy 
Hamiltonian is given as^ 

H(k) = ai |k|2 - a2(k ■ L)2 + je(B ■ L), (32) 

where L = (L^, Ly, L^) with 




/ 


0 0 Of 
0 0 1 
0-10, 


L,, = 



0 0 - 1 ) 
0 0 0 
10 0 , 


Lz = 



0 10 ) 
-10 0 
0 0 0 , 


Hu -H [-2fj cos( 7 r /4 -H 5) -i- 2t^ cos(; 7 r/ 4 )]crj, 
in other words. 


(28) This low-energy Hamiltonian harbors a pair of Weyl points 
at k = (0,0, ±ko), as we now show. Hereafter, we take B = 
(0,0, B) with pB > 0. By an explicit calculation, we have 



(FcTz 

0 

0 

0 ' 



( (a i- a 2 )kl 

-ipB 

0 ) 

Ho ^ Ho + 

0 

0 

Fo-z 

0 

0 

Fo-z 

0 

0 

(29) 

H((0,0,kz))^ 

ipB 

0 

(ai - a2)kl 
0 

0 

a\k\ , 


1 0 

0 

0 

Fo-z, 







(33) 


where F = - 2 fjCos( 7 r /4 -H 5) H- 2fj cos(; 7 r/ 4 ) With this 

term included, we can repeat the previous derivation of //eft 
and obtain the new low-energy Dirac Hamiltonian for ko - 
njA + 6: 

= (VxqxCTx + VyqyO-y + VzqzO-zTz) 

+{mcrzT+ + m*<JzT-) -I- (A -H Vz6)a-z^. (30) 


along the line k = (0,0,^^). The three eigenvectors are |1) = 
( 0 , 0 , 1 )^ | 2 ) = (l,/, 0 )^/V 2 , and |3) = ( 1 ,-/, 0 )^/V 2 , with 
energies 

E\ - aik^, 

E2 — (tti ~ Q'2)^j + PB, 

£3 = (ai - 02)14-PB. ( 34 ) 


The energy eigenvalues are 

£'i,2,3,4(q) = ± ^|vlql + v'lqj H- (A -I- Vz6 ± ^jv^q^ + \^)\3l) 


When k^ - ±ko = ± yjpBla 2 , we have the band touching 
El = E 2 , which defines the Weyl points. By tuning the param¬ 
eters of the gyroid photonic Weyl crystals, one can tune ko to 
(0, 0 , 7 r/ 4 ). Since £3 is far away from £i,£ 2 , we project the 
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Hamiltonian //(k) to the subspace spanned by |1) and |2). Let 
us introduce three Pauli matrices in this two-dimensional 
space. Near K± = (0,0, ±ko), we find that the effective two- 
band Hamiltonian reads 


//eft = ±kQ[—a2{qxO-x 

V2 


- q'yO-y) - a2qzCrz + (2a'i - a2)qz] 


+ ■■■(35) 


where qi = ki - (K±);, and “■ ■ ■” denotes higher order terms. 
Omitting the last term, which is proportional to the unity ma¬ 
trix, we have 

Hw ^Tz'^ ViCTiqi, (36) 


where Tz - refers to K±, Vx - —Vy - and = 

-a 2 ^o- Here, Ea. (l36l l takes almost the same form as Eq.(2) in 
the paper, except for the inessential difference that Tj appears 
in all the three terms (instead of appearing only in one of these 
three terms). The Dirac mass term for Ea.(l36ll is 

Hm - mT+ + m*T-, (37) 

where t± - (Tx ± iTy)l2. Compared to the Hm given in the 
paper, there is no cr^ factor here. This is due to the slight 
difference in the form of Hw- In each case, the Dirac mass 
term anticommutes with all terms in Hw- To generate such a 
Dirac mass, we only need to add the following term; 

H'^ ^mja-o(j), (38) 

j 


where ctq is the 2x2 unit matrix. Repeating the calculation in 
SecUof this Supplemental Material, we can obtain the Dirac 
mass 


m 


1 

^(mA 

1 


^{niA 


+ iniB — me — imo) 
+ + e^'^mc + 


(39) 


The simplest choice is (mA,mB,mc,mo) - 
(mo, mo, -mo, -mo), which leads to m = mo exp(/7r/4)/ V2. 

In the low-energy theory, the presence and absence of each 
term can be determined by symmetries. Erom symmetry con¬ 
sideration, any modulation of the photonic lattice with the A , B 
unit cells different from the C, D unit cells will produce the H' 
term in the low-energy effective theory, thus it will generate a 
Dirac mass, though the precise form of the band structure and 
the precise value of m can be determined only by photonic 
band structure calculations, which is beyond the scope of the 
present work. The existence of unidirectional modes along 
the vortex line of mass is nevertheless a topologically robust 
phenomenon. 

In conclusion, in this last section of the Supplemental Ma¬ 
terial, we have shown that the two-band model we have stud¬ 
ied in the paper (and its slight modification) captures the low- 
energy physics of the gyroid photonic Weyl materials, there¬ 
fore, our prediction of unidirectional modes along the vortex 
line of mass should be observable in these materials. 


